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Quantifying  absorption  by  aquatic  particles;  A  multiple  scattering  correction  for 
glass-fiber  filters 


ylbi/wc/— Absorption  spectra  measured  for  aquat¬ 
ic  particles  concentrated  onto  glass-fiber  filters  re¬ 
quire  a  correction  for  the  increase  in  pathlcngth  caused 
by  multiple  scattering  in  the  glass-fiber  filter.  A  mul¬ 
tiple  scattering  correction  was  calculated  from  op¬ 
tical  density  spectra  for  48  phytoplankton  cultures 
of  seven  species  representing  a  variety  of  cell  sizes, 
pigment  groups,  and  cell- wall  types.  The  relationship 
between  optical  density  in  suspensions  and  on  filters 
was  not  wavelength-dependent.  Differences  between 
blank  filters  were  always  spectrally  neutral.  Small 
differences  between  relationships  for  single  species 
were  inconclusive.  Given  the  absence  of  wavelength- 
dependent  effects,  we  report  a  single  general  qua¬ 
dratic  relationship,  OD,„,„(A)  =  0.378  OD6„(X)  + 
0.523  ODn„(X)^  (r^  =  0.988),  for  correcting  glass-fiber 
filter  spectra.  For  independent  samples,  the  average 
error  in  predicting  OD,„,„(X)  with  this  algorithm  at 
any  wavelength  was  2%.  Greatest  errors  were  in  spec¬ 
tral  regions  of  low  absorption.  Absorption  spectra 
for  particles  concentrated  onto  glass-fiber  filters  can 
be  quantitatively  corrected  for  multiple  scattering 
within  this  limit.  Applicability  of  the  algorithm  to 
field  samples  of  varied  composition  was  enhanced 
by  using  a  large  number  of  spectra  and  a  range  of 
cell  types  in  algorithm  development. 
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Variability  in  the  quantity  and  spectral  qual¬ 
ity  of  light  absorbed  by  particles  is  a  major  j 

cause  of  variability  in  light  attenuation  in  the  j 

ocean.  Unfortunately,  absorption  by  particles 
is  difficult  to  measure  because  suspended  par-  | 

tides  both  scatter  and  absorb  light  and  they  I 

are  present  in  low  concentration  in  much  of  , 

the  ocean.  Yentsch  (1962)  pioneered  a  tech-  j 

nique  that  increases  the  absorption  signal  to  a  , 

measurable  level;  he  concentrated  particles  i 

onto  a  membrane  filter  and  measured  their 
absorption  spectrum  while  the  particles  re-  j 
mained  on  the  filter.  Glass-fiber  filters,  which 
act  as  optical  diffusers  in  analogy  with  Shi-  ] 

bata’s  (1958)  opal  glass  technique,  are  now  ; 

commonly  used,  but  quantitative  use  of  these  j 

data  requires  correction  for  the  increase  in  ef¬ 
fective  pathlength  caused  by  multiple  scatter-  [ 
ing  within  the  filters.  This  technique  does  not  i 
measure  absorption  by  a  single  particle  but 
instead  represents  the  population  of  particles,  i  ^ 

Multiple  scattering  by  the  glass  fibers  in  the  -j - ^ 

filter  overshadows  the  increased  potential  for  V  i 

particle-particle  interactions  caused  by  con¬ 
centration  of  the  particles.  Subtleties  such  as 
variations  in  specific  absorption  coefficients 
that  arise  from  pigment  packaging  effects  are  -  -’-li 
still  apparent  in  spectra  measured  by  this 

method  (e.g.  Bricaud  and  Stramski  1990)  but  . 

these  variations  represent  the  population  av¬ 
erage.  Controversy  regarding  the  best  ap- 
proach  to  the  multiple  scattering  correction  re- 
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mains,  and  published  results  vary  (e.g.  Kishino 
et  al.  1985;  Maske  and  Haardt  1987;  Bricaud 
and  Stramski  1990;  Mitchell  1990).  This  note 
supports  the  use  of  an  empirical  quadratic  al¬ 
gorithm,  lacking  wavelength  dependence,  for 
this  correction.  Agreement  with  Mitchell’s 
(1990)  result  suggests  that  general  application 
of  a  single  algorithm  may  be  possible,  allowing 
use  of  consistent  methods  between  different 
investigators. 

Kiefer  and  SooHoo  (1982)  refined  the  filter 
method,  minimizing  loss  of  scattered  light  by 
placing  the  filter  close  to  the  detector  and  add¬ 
ing  a  correction  for  the  increase  in  effective 
pathlength  caused  by  scatter  within  the  glass- 
fiber  filter.  This  pathlength  amplification  fac¬ 
tor,  /3,  is  defined  as  the  ratio  of  the  optical  to 
geometrical  pathlength  (Butler  1962).  Kiefer 
and  SooHoo  (1982)  proposed  that  is  close 
to  6  for  two  layers  of  stacked  glass-fiber  filters 
(most  subsequent  investigators  use  a  single  lay¬ 
er).  After  correction  of  the  pathlength  with  0, 
absorption  coefficient  spectra  [ao(X);  m'‘]  are 
calculated  from  optical  density  spectra  (di¬ 
mensionless): 


2.3[OPn,.(X)  -  OPr, 1,(750)] 
^(V/A) 


(1) 


where  2.3  converts  from  log  base  1 0  to  log  base 
e,  V  (m’)  is  the  volume  filtered,  and  A  (m^)  is 
the  clearance  area  of  the  filter. 

Efforts  to  evaluate  0  have  increased  in  com¬ 
plexity  since  Kiefer  and  SooHoo’s  (1982)  as¬ 
sumption  of  constancy.  0  has  been  calculated 
from  comparisons  of  suspensions  and  filters 
and  shown  to  vary  with  the  optical  density  of 
the  sample  on  the  filter.  Mitchell  and  Kiefer 
(1988)  presented  a  nonlinear  relationship  be¬ 
tween  j9(X)  and  OPr,i,(X)  for  Whatman  GF/C 
glass-fiber  filters.  They  related  0{X)  for  other 
types  of  filters  to  the  equation  for  GF/C.  filters. 

Kishino  et  al.  (1985)  simply  calculated  the 
pathlength  amplification  factor  as  the  ratio  of 
filter  to  suspension  values,  without  considering 
0  to  vary  with  optical  density.  For  three  dif¬ 
ferent  cultures  of  phytoplankton  and  one  field 
sample,  0  (averaged  over  400-700  nm)  ranged 
from  2.43  to  4.7 1  (Kishino  et  al.  1985).  Kishi¬ 
no  et  al.  suggested  that  0  will  vary  with  particle 
density,  spectrophotometer  type,  and  filter 
type;  however,  their  presentation  of  data  av¬ 
eraged  over  the  visible  light  region  prevents 


analysis  of  ;3  as  a  function  of  optical  density 
in  the  manner  of  Mitchell  and  IGefer  (1988). 

Maske  and  Haardt  (1987)  discussed  path- 
length  amplification  caused  by  the  diffuse  na¬ 
ture  of  the  light  passing  through  the  sample 
and  filter  but  made  no  attempt  to  quantify  the 
pathlength  amplification  factor.  In  their  data, 
specific  absorption  coefficients  [absorption 
normalized  to  chlorophyll  concentration; 
ap*(675)]  decrease  as  chlorophyll  loading  on 
the  filter  increases;  increased  chlorophyll  load¬ 
ing  is  equivalent  to  an  increase  in  ODf,„(67  5), 
which  can  be  achieved  through  greater  volume 
filtered  or  greater  cell  density  in  the  same  vol¬ 
ume.  Because  average  values  of  0{X)  are  be¬ 
tween  2  and  2.5,  inclusion  of  a  pathlength  am¬ 
plification  factor  would  lower  their  ap*(675) 
by  this  factor.  Furthermore  since  the  nonlinear 
relationship  between  0iX)  and  ODm,(X)  is  steep¬ 
er  at  low  ODfii,(X)  (see  Mitchell  and  Kiefer 
1988),  application  of  0(X)  would  influence 
ap*(675)  at  low  filter  loads  more  strongly  than 
at  high  filter  loads.  For  low  ODf,„(X),  the  value 
of0(X)  will  be  higher  and  the  calculated  a  *(61 5) 
will  decrease  by  a  greater  amount  than  for  high 
ODf,„(X).  Application  of  a  nonlinear /3(X)  would 
reduce  the  observed  trend  in  <2p*(675)  with 
chlorophyll  loading. 

Yentsch  and  Phinney  (1989)  compared 
spectra  measured  for  suspensions  in  1 0-cm  cu¬ 
vettes  with  a  blank  glass-fiber  filter  acting  as 
an  optical  diffuser  to  spectra  obtained  on  fil¬ 
ters.  The  absorption  magnitude  measured  with 
these  two  techniques  was  the  same  (their  figure 
2)  so  they  did  not  apply  a  pathlength  ampli¬ 
fication  factor  to  their  absorption  data.  Meth¬ 
odological  differences  between  this  and  other 
approaches  may  underlie  the  contradictory  re¬ 
sults.  For  example,  a  10-cm  cuvette  with  an 
optical  diffuser  on  the  end  may  not  capture 
scattered  light  as  efficiently  as  an  integrating 
sphere  or  a  short  pathlength  cuvette  with  a 
diffuser.  Scattered  light  lost  along  the  sides  of 
the  10-cm  cuvette  may  never  reach  the  optical 
diffuser,  making  these  spectra  an  overestimate 
of  actual  absorption  and  explaining  the  ap¬ 
parent  absence  of  pathlength  amplification  in 
their  system. 

In  comparing  suspension  and  filter  values  of 
optical  density  for  three  algal  cultures,  Bricaud 
and  Stramski  (1990)  found  that  0(X)  varied 
significantly  at  low  ODf,|,(X)  and  observed  a 
“hysterisis  effect”  where  0(X)  was  not  the  same 
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for  identical  values  of  ODf,„(X)  within  a  par¬ 
ticular  spectrum.  For  ODfii,(X)  >0.2,  /3(X)  was 
less  variable.  They  corrected  their  field  data 
by  means  of  an  exponential  equation  for  /3(X) 
which  represented  an  average  relationship  de¬ 
veloped  from  their  laboratory  data. 

Using  a  Perkin-Elmer  Lambda  6  spectro¬ 
photometer,  Mitchell  (1990)  determined  that 
the  acceptance  angle  between  the  sample  and 
light  detector  did  not  affect  the  measured  op¬ 
tical  density,  i.e.  differences  in  optical  geom¬ 
etry  did  not  influence  the  relationship  between 
pathlength  amplification  and  optical  density. 
Mitchell  (1990)  assumed  no  wavelength  de¬ 
pendency,  produced  independent  nonlinear 
equations  for  four  filter  types,  and  advised  use 
of  a  direct  relationship  estimating  OD5u5p(X) 
from  ODfii,(X)  rather  than  calculation  of  /3(X) 
as  a  function  of  ODr,i,(X).  Absorption  coeffi¬ 
cient  spectra  are  calculated  by  subtracting 
ODf,i,(750)  from  OD,,|,(X),  estimating  OD,usp(X) 
with  the  quadratic  algorithm,  and  incorporat¬ 
ing  the  volume  filtered  and  the  filter  area; 


ap(X)  = 


2.30P,p,p(X) 

(V/A) 


(2) 


Mitchell  concluded  that  the  precision  of  the 
glass-fiber  filter  technique  is  10%,  the  accuracy 
is  1 5%,  and  variations  between  filter  lots  may 
affect  the  accuracy  of  this  method. 

Due  to  the  inconsistency  of  pathlength  am¬ 
plification  factors  reported  in  the  literature, 
skepticism  remains  regarding  the  validity  of 
correcting  absorption  spectra  measured  on 
glass-fiber  filters  for  multiple  scattering.  In  or¬ 
der  to  assure  the  proper  treatment  of  our  own 
data,  determine  whether  pathlength  amplifi¬ 
cation  for  our  spectrophotometer  agrees  with 
that  published  for  other  spectrophotometers, 
and  directly  develop  an  algorithm  for  the  GF/F 
filters  that  we  use  exclusively  in  our  field  and 
laboratory  sampling,  we  collected  data 'in  fall 
1989  for  determination  of  the  pathlength  am¬ 
plification  factor.  We  used  a  Kontron  Uvikon 
860  dual-beam  spectrophotometer  equipped 
with  a  scattered  light  transmission  accessory. 
This  accessory  positions  turbid  samples  close 
to- the  photomultiplier  tube  and  provides  a 
wide-acceptance  angle  for  capturing  scattered 
light. 

The  optical  density  of  algal  cultures  in  sus¬ 
pension  was  measured  in  1  -cm  cuvettes  with 
the  side  facing  the  light  detector  composed  of 


light-diffusing  opal  glass.  Because  of  their  small 
size  and  low  refractive  index  relative  to  water, 
algal  cells  primarily  scatter  in  the  forward  di¬ 
rection.  The  opal  ^ass  and  close  proximity  of 
the  detector  combine  to  capture  the  forward- 
scattered  light,  allowing  measurement  of  ab¬ 
sorption.  Bricaud  et  al.  (1983)  estimated  the 
loss  due  to  scattering  as  <0.5%  for  a  similar 
optical  orientation.  Optical  density  spectra  for 
suspended  particles  measured  in  these  opal 
glass  cuvettes  match  those  measured  with  an 
integrating  sphere  (Weidemann  and  Mitchell 
pers.  comm.).  Filtered  culture  medium  (Mil- 
lipore  GS,  0.22-Mm  pore  diameter)  was  used 
as  a  blank. 

The  optical  density  of  concentrated  cultures 
was  measured  after  filtration  of  a  subsample 
onto  a  Whatman  GF/F  glass-fiber  filter  (effec¬ 
tive  pore  size,  0.7  ^m)  at  low  vacuum.  Sub¬ 
sample  volume  was  chosen  so  that  the  geo¬ 
metrical  pathlength  of  filtered  samples  (volume 
filtered  divided  by  the  clearance  area  of  the 
filter)  matched  the  1-cm  pathlength  of  the  cu¬ 
vettes,  thus  ODfii,(X)  differed  from  (X) 

by  the  magnitude  of  the  pathlength  amplifi¬ 
cation  factor.  Prefiltered  seawater  (Millipore 
GS)  was  passed  through  a  second  GF/F  filter 
and  this  wet  filter  was  used  as  a  blank.  A  base¬ 
line  spectrum  for  two  wet,  blank  GF/F  filters 
was  stored  and  the  baseline-correction  feature 
of  the  spectrophotometer  used.  A  drop  of  fil¬ 
tered  seawater  adhered  the  wet  filters  to  a 
Plexiglas  sheet;  the  Plexiglas  was  then  placed 
in  a  custom-designed  holder  which  held  it  up¬ 
right  in  the  scattered  transmission  accessory. 
Thorough  saturation  of  both  filters  was  en¬ 
sured  by  letting  filtered  seawater  accumulate 
at  the  base  of  the  filter  holder;  however,  recent 
laboratory  measurements  indicate  that  optical 
density  and  spectral  characteristics  of  the  sam¬ 
ple  do  not  vary  with  water  content  (Cleveland 
unpubl.).  The  sample  side  of  the  filter  faced 
the  light  source  and  the  Plexiglas  faced  the 
detector.  OD,u5p(X)  and  ODfii,(X)  were  mea¬ 
sured  from  760  to  390  nm  at  4-nm  spectral 
bandwidth. 

A  comparison  of  blank  GF/F  filters  within 
the  same  lot  (n  =  6)  and  between  lots  («  =  7) 
was  made.  The  baseline  was  stored  with  two 
filters  from  the  same  lot,  and  each  successive 
filter  was  compared  to  the  same  reference  filter. 
ODr,i,(X)  for  all  “sample”  blank  filters  was  spec¬ 
trally  flat  with  respect  to  the  reference  filter. 
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Four  of  these  spectra  showed  a  slight  elevation 
of  ODfi„(400)  (0.004  relative  optical  density 
units)  compared  to  ODni,(750),  with  a  smooth, 
gradual  slope  between  400  and  750  nm.  For  a 
particulate  sample  with  optical  density  of  in¬ 
termediate  value  (e.g.  0.2  optical  density  units), 
this  amounts  to  a  2%  change  across  the  spec¬ 
trum.  This  variability  is  insignificant.  Spectra 
were  always  flat  but  optical  density  was  not 
always  zero.  However,  this  difference  is  ac¬ 
counted  for  when  ODf.n(750)  is  subtracted  from 
ODfi,.(X). 

Unialgal  phytoplankton  cultures  were  grown 
in  sterile  filtered,  nutrient-enhanced,  artificial 
seawater  at  20®C  under  “cool-white”  fluores¬ 
cent  lamps  (33  ^mol  photons  m“^  s  ')  on  a  12 : 
12  L/D  cycle.  Cultured  phytoplankton  were 
used  because  natural  samples  are  not  dense 
enough  to  allow  direct  measurement  of  optical 
density  on  unconcentrated  samples.  Algal  spe¬ 
cies  ranged  from  about  1  to  20  /im  in  diameter. 
We  measured  48  pairs  of  spectra  for  seven 
different  species:  Chaetoceros  gracilis  ( 1 2  pairs 
of  spectra),  Thalassiosira  weisspogii  (11),  Am- 
phidinium  carterae  (5),  Micromonas  pusitta  (2), 
Prasinophyte  fi  48-23  (3),  Dunalielta  tertiolec- 
ta  (1 1),  and  Synechococcus  WH7803  (4).  Each 
species  was  grown  to  a  variety  of  cell  densities, 
providing  a  range  of  optical  densities.  For  rep¬ 
licate  samples,  the  C.V.  ranged  from  1  to  5% 
for  suspensions  and  1  to  9%  for  fitters.  Highest 
variation  occurred  at  low  optical  density.  Spec¬ 
tra  were  measured  immediately  after  filtration. 
No  qualitative  degradation  effects,  such  as  the 
blue-shifted  Soret  peaks  described  by  Stramski 
( 1 990),  were  observed  and  precautions  against 
pigment  degradation  were  not  taken.  As  point¬ 
ed  out  by  Stramski,  only  certain  species  exhibit 
degradation  effects,  and  our  species  are  not 
among  these.  He  recommends  a  systematic 
study  of  the  effect  of  paraformaldehyde  rather 
than  general  adoption  of  presetwatives  as  a 
routine  technique. 

The  respective  optical  densities  at  750  nm 
were  subtracted  from  OD,u,p(X)  and  OD„„(X). 
Spectra  measured  on  glass-fiber  filters  are  fea¬ 
tureless  in  this  region,  despite  the  temperature 
dependent  change  that  occurs  in  water  ab¬ 
sorption  at  745  nm  (Pegau  and  Zaneveld  1 993); 
filters  apparently  do  not  contain  enough  water 
to  exhibit  this  effect.  Examination  of  spectra 
for  all  48  suspensions  measured  in  1-cm  cu¬ 
vettes  also  did  not  reveal  any  features  near  745 


nm.  We  have  observed  temperature-depen¬ 
dent  artifacts  with  10-cm  cuvettes  in  the  Kon- 
tron  spectrophotometer  (see  Pegau  and  Za¬ 
neveld  1993)  as  the  reference  cell  warms  in 
the  sample  compartment.  Similar  artifacts  did 
not  occur  in  the  1  -cm  cuvettes  used  here.  The 
same  techniques  were  used  with  both  1 0-  and 
1-cm  cuvettes,  indicating  that  the  1-cm  cu¬ 
vettes  had  sufficient  time  in  the  sample  com¬ 
partment  to  achieve  a  temperature  differential, 
but  the  signal  was  not  large  enough  to  be  mea¬ 
sured. 

Data  analysis  was  restricted  to  the  visible 
light  region,  400-700  nm.  Below  400  nm, 
OD,,„(\)  is  noisy  due  to  low  light  through-put. 
Values  of  ODni,(X)  >0.4  were  eliminated  from 
the  data  set  because  our  field  and  laboratory 
data  are  always  lower  than  0.4  and  generally 
below  0.3.  Restriction  of  the  data  range  also 
provides  consistency  with  the  procedure  used 
by  Mitchell  (1990). 

Qualitatively,  results  for  /3(\)  as  a  function 
of  ODf,i,(X)  (not  shown)  resembled  those  of  pre¬ 
vious  studies.  The  mean  value  of  /3(X)  was  2.46, 
/8(X)  varied  as  an  inverse  function  of  ODf,|,(\), 
and  |3(X)  was  more  variable  at  low  ODfi|,(X). 
The  “hysteresis  effect”  observed  by  others  was 
not  apparent  in  all  pairs  of  spectra  and  was 
small  when  present,  as  opposed  to  the  large 
loops  shown  by  Bricaud  and  Stramski  (1990). 

Instead  of  analyzing  0(X)  as  a  function  of 
OD,,„(X).  we  followed  the  suggestion  of  Mitch¬ 
ell  (1990)  and  examined  OD5u,p(X)  as  a  func¬ 
tion  ofODr,i,(X)  (n  =  13,256).  Statistical  anal¬ 
yses  were  performed  with  International 
Mathematics  and  Statistical  Libraries  (IMSL). 

goodness-of-fit  tests  indicated  that  a  second- 
order  polynomial  represented  the  data  better 
than  a  third-order  polynomial  or  straight  line. 
The  resulting  quadratic  equation, 

OD,„,p(X)  =  0.378  ODfi„(X) 

+  0.523  ODn,.(X)^  (3) 

(r^  =  0.988),  calculates  OD,u,„(X)  from  mea¬ 
sured  ODn„(X).  The  95%  confidence  intervals 
are  0.371-0.384  for  the  first  coefficient  and 
0.496-0.550  for  the  second.  An  intercept  was 
not  included  because  absorption  by  the  sus¬ 
pension  should  equal  zero  when  absorption  by 
the  filtered  sample  equals  zero.  This  algorithm 
is  statistically  valid  only  for  the  data  range 
[ODr,ii(X)  between  0  and  0.4]  over  which  it  was 
developed. 
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Values  of  ODsu5p(X)  predicted  from  Eq.  3 
differ  by  <  10%  from  the  estimates  of  Mitchell 
( 1 900),  with  the  curves  diverging  to  this  max¬ 
imum  as  ODfinfX)  increases.  This  result  is  with¬ 
in  the  expected  accuracy  of  the  method,  even 
though  coefficients  in  each  algorithm  do  not 
fall  within  the  95%  confidence  intervals  for  the 
other.  Differences  in  optical  geometry  for  var¬ 
ious  spectrophotometers  may  influence  the  re¬ 
lationship  between  ODsp5p(X)  and  ODnufX). 
Mitchell  suggested  that  filter  thickness  varies 
with  filter  lot,  which  would  also  affect  the  re¬ 
lationship  between  ODspjpfX)  and  ODr,i,(X).  Our 
data  were  obtained  before  his  were  published; 
various  filter  lots  were  used  and  lot  number 
was  not  recorded.  However,  as  described  pre¬ 
viously,  GF/F  blanks  from  seven  different  lots 
showed  no  spectral  differences,  and  differences 
in  magnitude  between  various  lots  have  no 
effect  on  the  algorithm  or  the  sample  absorp¬ 
tion  after  subtraction  of  ODr,i,(750).  Despite 
differences  in  coefficients,  the  two  algorithms 
give  similar  results,  suggesting  that  use  of  a 
single  algorithm  for  various  spectrophotome¬ 
ters  may  be  possible. 

Accuracy  of  the  algorithm  was  examined 
with  10  independent  samples  with  ODr,i,(435) 
between  0.03  and  0.3.  OD,p5p(X)  predicted  from 
ODr,i,(X)  with  Eq.  3  was  within  2%  of  directly 
measured  OD,p,p(X)  when  averaged  over  all 
wavelengths  and  samples,  with  a  maximum 
error  of  20%  for  one  sample  in  the  yellow 
wavelength  region  where  absorption  is  low.  The 
optical  densities  of  the  filtered  test  samples 
were  restricted  to  the  range  (0-0.4)  used  for 
algorithm  development  and  typical  of  field 
samples  by  controlling  the  volume  filtered. 

The  potential  for  wavelength  dependence  of 
the  multiple  scattering  correction  was  exam¬ 
ined  by  fitting  a  second-order  polynomial  at 
30  separate  wavelengths  [every  10  nm,  n  =  48 
except  at  440  nm  where  two  values  of 
ODf, 1,(440)  exceeded  0.4].  Equations  for  wave¬ 
lengths  from  400  to  480  nm  were  identical. 
Differences  between  this  group  and  other 
wavelengths  did  not  follow  any  pattern  with 
wavelength,  as  would  be  expected  from  phe¬ 
nomena  such  as  wavelength-dependent  scat¬ 
tering.  Instead,  these  differences  appeared  to 
be  an  artifact  of  inherent  differences  in  the 
range  of  optical  density  at  each  wavelength. 
Phytoplankton  absorption  spectra  have  peaks 
and  valleys  at  consistent  wavelengths,  causing 


the  maximum  value  of  ODni,(X)  to  reach  0.4 
at  440  nm  but  only  0.09  at  560  nhi.  As  a  con¬ 
sequence,  the  second  coefficient,  C2(X),  of  the 
wavelength-specific  second-order  polynomials 
varied  inversely  with  mean  ODr,u(X)  at  the  re¬ 
spective  wavelengths  (H  =  0.71).  The  first  co¬ 
efficient  of  the  polynomials,  C|(X),  varied  less 
and  did  not  show  a  similar  trend  with  range 
of  ODf,„(X).  C,(X)  for  only  5  of  the  30  wave¬ 
lengths  did  not  fall  within  the  95%  C.I.  for  all 
other  wavelengths,  indicating  that  spectral  dif¬ 
ferences  were  minimal. 

Wavelengths  with  similar  ranges  of  ODf,i,(X) 
(e.g.  400,  480,  and  680  nm.  Fig.  1)  had  statis¬ 
tically  indistinguishable  algorithms  (t-test,  P 
<  0.05)  even  when  separated  by  many  nano¬ 
meters,  such  as  the  280-nm  difference  between 
400  and  680  nm.  Because  wavelength-specific 
algorithms  were  similar  when  ranges  of  ODf,i,(X) 
were  similar,  we  ascribe  the  apparent  wave¬ 
length  dependence  to  a  statistical  artifact  caused 
by  inherent  differences  in  the  range  of  data  at 
the  various  wavelengths.  No  trends  in  the  error 
between  predicted  and  measured  OD^^spfX)  as 
a  function  of  wavelength  were  apparent,  i.e. 
no  consistent  patterns  of  over-  or  underesti¬ 
mation  at  a  particular  wavelength  occurred, 
verifying  the  absence  of  wavelength-depen¬ 
dent  trends  in  the  algorithm. 

Spectral  differences  attributable  to  scattering 
by  glass-fiber  filters,  as  described  by  Roesler 
and  Perry  (unpubl.),  were  not  evident  in  these 
data.  If  present,  spectral  scattering  by  glass- 
fiber  filters  would  appear  as  increases  in  ap¬ 
parent  absorption  at  the  blue  end  of  the  spectra 
or  as  increases  in  the  ratio  of  ODf, ,,(400): 
ODf,„(675)  compared  to  OD5ujp(400): 
OD,u,p(675).  Even  for  spectra  with  low  optical 
density,  where  these  changes  would  be  more 
obvious,  no  increase  or  flattening  of  ODfii,(400) 
relative  to  OD5u,p(400)  was  present.  Ratios  of 
ODf, I, (400) :  ODfii,(675)  for  all  48  pairs  of  spec¬ 
tra  used  in  algorithm  development  were  evenly 
divided  between  greater  and  lower  than  the 
corresponding  ratio  for  the  suspension.  We 
speculate  that  blank  fillers  more  closely 
matched  sample  filters  in  our  study  compared 
to  Roesler  and  Perry’s  (unpubl.)  data.  Com¬ 
parison  of  GF/F  fillers  wiihin  and  beiween  lols 
showed  no  speclral  fealures,  again  discounting 
differences  in  spectral  scattering  between  fil¬ 
ters. 

Species-specific  algorithms  were  compared 
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Fig.  1.  Wavelength-specific  data  and  curve  fits  for  OD,„p(X)  as  a  function  of  OD|i„(\)  at  400  nm  (+,  solid  line), 
480  nm  (O,  dashed  line),  and  680  nm  (x,  dotted  line). 


for  those  species  with  sufficient  number  of 
spectra  and  range  of  optical  density  to  provide 
statistically  useful  algorithms.  Curves  for  the 
two  diatoms  with  siliceous  frustulcs,  C  gracilis 
(~6  ^m  equivalent  spherical  diameter)  and  T. 
weissjlogii  (~  1 1  nm),  were  similar.  A.  carterae 
(an  unarmored  dinoflagellate;  ~20  Mitt)  had  a 
steeper  relationship  (higher  Cj),  while  D.  ter- 
tiolecta  (a  naked  chlorophyte;  ~8  ttni)  had  a 
shallower  relationship  (lower  C,  and  C^.  These 
differences  are  intriguing  but  the  data  set  is  not 
diverse  enough  to  establish  trends  caused  by 
cell-wall  composition  or  cell  size.  The  species 
were  chosen  to  produce  a  general  algorithm 
rather  than  provide  a  close  examination  of  spe¬ 
cies  differences. 

The  analysis  presented  here  indicated  that 
the  correction  for  multiple  scattering  in  glass- 
fiber  filters  was  independent  of  wavelength,  al¬ 
though  species  independence  was  not  conclu¬ 
sively  shown.  Development  of  the  algorithm 


from  a  data  set  including  species  with  varied 
characteristics  and  a  large  number  of  data 
points  increases  the  generality  of  the  algorithm 
for  use  with  diverse  field  samples.  Tests  on 
independent  samples  showed  good  predictive 
ability.  Optical  density  spectra  measured  on 
glass-fiber  filters  can  be  quantitatively  cor¬ 
rected  with  this  algorithm. 
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